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Summary:  The  purpose  of  the  lecture  is  to  demonstrate  the  method  of  Matched  Asymptotic  Expansion  (MAE) 
to  be  efficiently  applicable  to  various  supercavitating  flow  problems.  Two  model  problems  were  chosen  to  prove 
the  statement.  The  first  one  is  a problem  of  the  flow  around  a supercavitating  shock  free  hydrofoil  with  spoiler 
mounted  on  its  trailing  edge  and  wedge-like  fully  wetted  leading  edge.  A gallery  of  local  nonlinear  edge  solution 
is  proposed  for  the  inner  flow  problems  in  the  vicinity  of  the  spoiler  as  well  as  in  the  vicinity  of  the  sharp  or 
rounded  cavitating  leading  edge.  The  second  problem  is  supercavitating  wing  of  large  aspect  ratio  beneath 
the  free  surface.  In  both  cases  analytical  solutions  to  nonlinear  inner  problem  were  applied  to  accomplish  the 
matching  procedure  and  to  significantly  extend  the  range  of  the  flow  parameters  where  asymptotic  approach 
gives  reliable  results.  Solutions  obtained  in  the  framework  of  the  MAE  method  as  well  as  linear  and  nonlinear 
theories  are  illustrated  in  the  text  by  numerical  data. 

1.  Introduction 

The  Matched  Asymptotic  Expansions  (MAE)  method  is  an  appropriate  solution  technique  for  a wide  variety 
of  fluid  mechanics  problems.  As  a rule,  the  classical  MAE  method  involves  the  following  steps  while  solving  a 
lifting  flow  (including  cavitating  one)  problem: 

a)  to  determine  what  the  small  parameters  of  the  problem  are  and,  depending  on  the  answer  to  the  question, 
subdivide  the  flow  domain  into  a so-called  farfield  and  a.  nearfield,  that  is,  regions  far  from  and  in  the  vicinity 
of  a source  of  singular  perturbations;  formulation  of  problems  in  those  regions  - outer  and  inner  problems 
correspondingly.  The  main  goal  of  the  subdividing  procedure  is  to  simplify  the  outer  problem  as  compared  to 
a general  one  (for  instance,  such  a simplification  is  brought  about  by  the  linearization  procedure  for  all  the 
boundary  conditions  in  outer  region),  and  to  take  into  account  as  much  nonlinear  effects  as  possible  in  the 
nearfield. 

b)  a ‘rough’  asymptotic  analysis  of  both  outer  and  inner  descriptions  from  the  viewpoint  of  a solution  class 
to  the  problems  arising  in  those  subdomains.  For  instance,  a presence  of  the  spoiler  on  the  trailing  edge  of 
planing  or  cavitating  hydrofoil  dictates  a non-traditional  class  (namely,  oo  — oo  instead  of  customary  oo  — 0) 
of  linearized  outer  solution  to  the  corresponding  mixed  boundary  problem.  In  the  case  of  a cavitating  problem 
that  means  that  the  function  of  complex  conjugate  velocity  in  the  farfield  has  to  have  square  root  singularity 
on  the  trailing  edge  of  the  hydrofoil  and  —1/4  on  the  leading  one,  see  Rozhdestvensky  & Fridman  [17]; 

c)  solution  of  the  outer  lifting  flow  problem,  in  the  flow  region  far  from  the  inner  zone.  Such  an  outer  asymptotic 
expansion  is  often  derived  under  the  assumption  of  small  perturbations  brought  by  a hydrofoil  into  the  inflow. 
The  assumption  loses  its  correctness  in  the  vicinity  of  the  inner  sub-region  and  so  does  the  outer  expansion. 
Outer  solution  contains  some  elements  of  uncertainty  caused  by  the  influence  of  local  problem,  as  some  significant 
features  were  omitted  from  consideration. 

Note  that  the  MAE  method  does  not  always  require  the  outer  problem  to  be  linear  and  it  is  not  a general 
restriction  on  its  applicability.  Examples  as  such  are  abound  but  here  we  to  some  degree  confine  ourselves  by 
such  an  assumption. 

d)  solution  of  the  inner  flow  problem,  i.e.  construction  of  a so-called  inner  asymptotic  expansion  which  is 
correct  in  the  local  region  and  loses  its  correctness  in  the  farfield.  It  should  be  pointed  out  that  inner  problem 
is  usually  considered  in  stretched  local  coordinates,  the  stretching  factor  being  connected  with  a chosen  small 
parameter.  It  is  by  taking  advantage  of  the  stretching  procedure  that  the  inner  problem  is  also  simplified,  for 
instance,  reducing  of  three-dimensional  problem  to  two  dimensions  in  the  nearfield,  neglecting  the  gravity,  etc. 
This  allows  to  apply  a nonlinear  approach  to  adequately  describe  the  most  important  part  of  the  flow  domain. 
Because  of  the  influence  of  the  outer  region,  inner  solution  contains  some  unknown  parameters,  as  well; 

e)  matching  procedure  for  outer  and  inner  expansions  which  is  applied  to  ‘blend’  them  and  to  take  into  account 
their  interaction.  This  stage  allows  us  to  close  the  solution  to  the  whole  problem  because  all  the  unknown 
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parameters  in  the  far-  and  nearfield  are  determined.  Using  the  results  obtained,  a composite  uniformly  valid 
asymptotic  solution  to  the  problem  under  consideration  is  reached.  Such  a composite  solution  is  shown  to 
be  correct  in  the  whole  flow  domain  and  enables  one  to  remove  some  intrinsic  drawbacks  of  the  purely  linear 
description  used  in  the  outer  region. 

The  advantages  of  such  an  asymptotic  method  of  treating  the  lifting  flow  problems  are  easily  seen  and  thor- 
oughly discussed  by  many  authors,  see  Van  Dyke  [24],  Cole  [2],  Ogilvie  [14],  Rozhdestvensky  [16],  etc.  In  this 
connection  a ‘mathematical  constructor’  approach  should  be  mentioned,  advocated  in  Fridman  [3]  and  Fridman 
& Rozhdestvensky  [4,  17].  It  is  by  applying  this  approach  that  one  can  pick  up  and  then  to  combine  the 
necessary  outer  and  inner  asymptotic  expansions  (just  like  standard  ‘mathematical  parts’  or  the  details  of  the 
whole  construction)  to  obtain  a desired  everywhere  correct  solution. 

It  is  a well  known  fact  that  the  MAE  method,  being  applied  to  the  lifting  flow  problems,  enables  one  to 
scrutinize  those  flow  regions,  where  nonlinearities  of  the  problem  under  consideration  are  concentrated  and 
most  pronounced,  while  outer  flow  domain,  far  from  such  zones  of  nonlinearity,  can  be  described  more  or  less 
sufficiently  in  the  frameworks  of  an  appropriate  linear  theory.  In  this  case  the  outer  description  is  certainly 
to  retain  some  aftermaths  of  the  linearization  of  the  boundary  conditions  which,  though,  can  be  efficiently 
overcome  by  MAE  technique. 

2.  Super cavitating  shock  free  hydrofoil  with  spoiler 

2.1  Problem  formulation 

Consider  a cavitating  problem  for  the  hydrofoil  with  the  spoiler,  at  arbitrary  cavitation  number,  see  figure  1,  the 
influence  of  gravity  being  neglected.  The  Efros  cavity  closure  model  with  re-entrant  jet  is  adopted.  The  origin 
of  the  Cartesian  coordinate  system  is  taken  at  the  plate’s  leading  edge,  a:-axis  being  directed  downstream  and 
y upwards.  There  is  an  incident  stream  with  speed  Voo  coming  from  the  left.  The  region  occupied  by  the  fluid 
is  bounded  by  the  solid  straight  boundaries  [AC],  [OB]  and  arbitrarily  arched  portion  (CO)  given  in  the  form 
y = f(x)  and  the  cavity  surfaces  (AE)  and  (BE),  the  intervals  [AC]  and  [OB]  being  of  the  length  [AC \ = lw 
and  | OB | = e.  The  hydrofoil  chord  |OC|  is  chosen  to  be  l.  The  incidence  angle  is  a,  the  inclination  angle  of 
the  spoiler  with  respect  to  the  hydrofoil  chord  is  (i  and  the  angle  made  by  [AC]  and  (CO)  is  7. 

It  is  of  importance  that  the  dividing  streamline  would  meet  the  cavitating  plate  at  the  vertex  of  the  wetted 
surface  (AC OB),  namely  at  point  C,  to  provide  the  shock-free  cavitating  mode.  In  this  case  the  length  of 
the  upper  wetted  portion  of  the  leading  edge  Zw  is  to  be  treated  as  an  unknown  parameter  of  the  problem,  the 
angle  7 made  by  [AC]  and  (CO)  being  given.  Since  the  cavitation  number  a > 0,  the  velocity  absolute  value 
011  the  free  surfaces  is  vo  = Voo  >/l  + a . The  point  at  infinity  on  the  re-entrant  jet  is  denoted  by  E and  is  quite 
distinct  from  D,  the  point  at  infinity  in  the  main  stream.  An  additional  stagnation  point  F appears  in  the  flow 
pattern  because  of  the  re-entrant  jet  influence.  That  is  why  a double-sheeted  Riemann  surface  is  introduced, 
one  sheet  carrying  the  main  flow  pattern  and  the  second  including  the  jet.  The  direction  of  the  jet  at  point  E 
is  determined  by  the  angle  fi , see  figure  1 . All  the  flow  parameters  are  rendered  nondimensional  by  using  chord 
length  l and  velocity  Voo. 


Figure  1:  Shock  free  supercavitating  hydrofoil  with  spoiler. 

As  is  customary,  the  problem  is  considered  to  be  solved  when  the  velocity  potential  function  ip(x,  y)  is  found.  The 
harmonic  function  <p(x,y),  being  the  real  part  of  an  analytical  function  of  complex  potential  w = ip  + iip,  has  to 
satisfy  the  boundary  kinematic  (flow  tangency)  condition  on  the  wetted  length,  the  dynamic  (pressure/velocity 
constancy)  condition  on  the  free  surface  and  the  condition  at  infinity. 
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Under  assumption  that  the  hydrofoil  brings  small  perturbations  into  the  inflow,  subdivide  the  flow  domain  into 
a farfield  (at  distances  of  the  order  of  0(1))  and  nearfields  (in  the  vicinity  of  the  edges)  we  apply  the  MAE 
method. 


2.2.  Outer  linear  solution 

Assume  that 


6(x)  = f'{x)  = o(l)  , a = o(l)  , e = y=o(l),  y = o(a),  <7  = o(l), 

that  is,  the  hydrofoil  and  the  cavity  brings  small  perturbations  into  the  inflow  and  therefore  the  linearization 
procedure  has  to  be  accomplished  for  the  nonlinear  boundary  (dynamic  and  kinematic)  conditions  both  on 
the  free  surfaces  and  on  the  wetted  portion  of  the  foil.  Under  those  circumstances  we  can  neglect  the  second 
order  terms  in  all  the  equations  and  can  formulate  the  linearized  flow  problem  for  supercavitating  contour  with 
spoiler  [22,  21,  9,  10,  17].  All  the  wetted  surface  of  the  hydrofoil  and  the  cavity  appears  as  a slit  of  length  L 
in  the  linearized  plane  z = x + i y,  where  L actually  is  a linearized  cavity  extent.  Boundary  conditions  for  real 
and  imaginary  parts  of  the  conjugate  velocity  \°  are  projected  on  the  upper  and  lower  boundaries  of  the  slit, 
see  figure  2,  where 

X°(z)  = u°  — iv°  = , 

v0  dz 

and  w(z)  = tp  + i'il>  - function  of  the  complex  potential. 


L — z 


1 


Figure  2:  Physical  and  auxiliary  planes. 


The  physical  2-plane  has  to  be  mapped  into  the  auxiliary  upper  semi-plane  £ where  the  mixed  boundary  value 
problem  is  formulated  in  oo  — oo  class  due  to  the  presence  of  the  spoiler  located  at  the  trailing  edge  of  the  foil. 
The  class  oo  — oc  means  that  the  function  x°(z)  has  to  have  square-root  singularities  for  2 = 0 and  2=1.  Such 
a behaviour  of  the  outer  solution  dictates  by  the  asymptotic  analysis  of  the  flow  problems  arising  in  the  vicinity 
of  the  leading  and  trailing  edge  with  the  spoiler  done  in  the  section  that  follows.  The  analysis  also  reveals  that 
lw/l  < a and  therefore  the  kinematic  boundary  condition  on  the  upper  wetted  portion  [AC]  of  the  leading 
edge  can  be  neglected  correctly  in  the  framework  of  the  linear  theory  for  all  the  wetted  surface  of  the  hydrofoil 
appears  as  a segment  of  unit  length  on  the  lower  boundary  of  the  slit  in  the  physical  2-plane  (points  A and  C 
coincide  and  so  do  points  B and  O,  and  E and  F). 

Note  that  in  the  case  of  lw/l  = 0(1),  the  kinematic  boundary  condition  on  the  upper  wetted  portion  [AC] 
has  to  be  satisfied.  To  ensure  shock  free  regime  of  the  flow  at  the  leading  edge,  the  solution  to  the  problem 
should  behave  at  point  C like  *Jz  as  2 — > 0.  Such  a solution  exists  for  unique  value  of  the  parameter  lw/l  which 
becomes  an  unknown  of  the  linear  problem. 

Back  to  the  problem  as  lw/l  <C  a,  apply  Keldysh-Sedov  formula  for  the  upper  semi-plane  in  oo  — oo  class: 


x°(C) 


iB 


v/CCChh 


+ i 


£ + a / 1 


a 

I 


t,  9(t)  dt 
a — t C + t 


+ A £ + C 


(1) 
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where  A,  B.C  and  a = 1 /y/L  - 1 are  unknown  parameters  to  be  determined  from  the  cavity  closure  condition 
and  condition  at  infinity 

Im  ^ X°(C)  = 0 ; X°(i)  = l-f, 

<=< 

along  with  condition  provided  by  the  matching  procedure  in  the  vicinity  of  the  spoiler.  It  is  to  be  underlined 
that  the  solution  provides  a square-root  singularity  for  y°  at  cavity  trailing  edge  z = L,  £ — > oo.  It  corresponds 
to  the  linearized  closed  cavity  closure  models:  Riaboucliinsky,  Efros-Kreisel-  Gilbarg,  Tulin-Terentev  (single 
spiral  vortex). 

Note  that  the  simple  non-quadrature  approach  can  be  proposed  [21]  for  a wide  range  of  the  functions  y = /( x) 
characterizing  the  lower  surface  distribution  of  the  hydrofoil. 

It  is  obvious  that  the  new  function 

where  9(x)  = f'(x)  is  a tangential  angle  to  the  foil  at  point  x,  has  pure  real  values  on  the  wetted  portion  (as 
x 6 [0;  1 ],  y = 0+)  and  pure  imaginary  values  on  the  cavity  surfaces,  see  figure  2.  Let  us  assume  that  the  wetted 
portion  of  the  cavitating  hydrofoil  is  a polynomial 


n 

f 0)  = J2aixi 

2=0 


and  therefore 

n 

9(x)  = f(x)  = '£iaixi~1 . 

2 = 1 

Then  function 

n 

fi(C)  = .\°(C)-l  + i X>.V_1(C) 

2 = 1 

where  z(£)  = LC,2/{  1 + £2),  has  to  satisfy  homogeneous  boundary  conditions  on  the  upper  semi-plane  £ in 
oo  — oo  class.  It  should  be  pointed  out  that  the  multiplicity  of  a pole  at  the  infinity  2 — > oo  for  the  function 
Q(z)  is  equal  to  (n  — 1),  i.e.  to  the  degree  of  the  polynomial  f{x).  That  is  why  the  solution  can  be  derived 
without  an  integration  procedure  and  is  of  the  form 

*°(c)=1+^wbriW1/¥ 


AC  + C + ^2 


T>k  + i£k  T^k 


i£k 


(C  -i)k  (C  +i)k 


(2) 


Condition  at  infinity  z = oo  and  at  its  image  £ = i for  function  x°(C) 

iim(x°(C)  - l)  = ; 


lim{  (x°(C)  - 1)  (C  - i)"'}  = 0 for  k = 1, 


, n - 1 . 


linearized  cavity  closure  condition  re-written  in  the  form 


dy° 


along  with  the  matching  condition  allow  unknown  parameters  A.  B,  C,  T>k-  £k  and  cavity  extent  L in  (2)  to  be 
determined.  So,  the  number  2 n + 2 of  unknowns  coincides  with  that  of  conditions. 

Once  function  y°(£)  is  found  (i.e.  the  unknowns  of  the  outer  problem  are  derived),  the  cavity  volume  is  deter- 
mined 

X X 

ZC+avO)  = j dt  ’ ^cavO)  = I"  + 1 “ i/f1)  > (3) 

0 1 

where  superscripts  + and  — denote  the  upper  and  lower  boundary  of  the  cavity  correspondingly. 


19-5 


All  the  hydrodynamic  coefficients  are  also  connected  to  x°(z): 


c-w--2 


(4) 


Cl  = I Cp(x)dx  = 2 Re  j{x°{ Q - 1)  ^ dC  , 

0 i 


CD  = 


(5) 

(6) 


Both  classical  and  nonquadrature  solutions  to  the  supercavitating  problem  remain  valid  for  the  open  cavity 
closure  scheme  (linearized  analog  of  the  Wu-Fabula  model)  as  well.  Such  a model  specifies  the  velocity  field 
be  continuous  at  the  trailing  edge  of  the  cavity  what  means  that  the  cavity  and  the  trailing  wake  conjugates 
smoothly.  From  the  mathematical  point  of  view  it  corresponds  to  non-singular  behaviour  of  x°(z)  at  point 
z = L which  implies  the  condition  A = 0 to  be  satisfied.  Then  the  number  of  unknowns  becomes  2n  + 1 and 

dy° 

the  cavity  closure  condition  Im  — j—  (i)  = 0 should  be  neglected. 


2.3.  Inner  nonlinear  descriptions 

2.3.1.  A simplest  spoiler  problem.  It  was  already  mentioned  above  that  outer  asymptotic  expansion  loses 
its  correctness  in  the  vicinity  of  the  spoiler  as  z — ► 1 and  inner  description  dictates  the  solution  class  of  the  outer 
one,  namely  oo  — oo.  The  picture  of  the  flow  shown  in  figure  3 corresponds  to  stretching  of  local  coordinates  in 
this  region  by  a factor  1/e:  X = (x  — l)/e,  Y = y/e,  e — > 0.  The  region  occupied  by  the  fluid  is  bounded  by  the 
solid  boundaries  [OB)  and  [OA]  and  the  free  surface  ( AB ),  the  interval  [OA]  being  of  the  unit  length  \OA\  = 1. 
The  absolute  value  of  velocity  at  ‘local’  infinity  and  on  the  free  surface  (AB)  is  an  unknown  parameter  to 
be  determined  through  the  matching  procedure  with  outer  problem.  This  is  the  simplest  flow  problem  for  a 
straight  spoiler. 


Figure  3:  Flow  pattern  and  auxiliary  /-plane  for  the  inner  spoiler  problem. 


With  the  correspondence  between  physical  and  auxiliary  planes  (the  first  quadrant  of  the  / = / + ijj  plane  is 
chosen  as  an  image  of  the  region  occupied  by  the  fluid),  one  can  easily  obtain  the  analytical  inner  nonlinear 
solution  (following  to  Chaplygin  method  of  singular  points)  in  the  form: 


(s) 


dw 

vJJz 


= N /; 


dw 

dl 


(7) 
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wherefrom 


*<•>=£( 


iv  /•<  /c  + 1 
s'c^i 


/3/tt 


dC, 


where  iV  - unknown  parameter  to  be  determined  from  the  condition  Za  = Z( 0)  = e 1/3  which  yields 


N 

vL 


1 /l+£\£'/7r 


IX 


i-i 


The  pressure  distribution  coefficient  (0  < £ < oo)  is 

O’-i  - |xi(C)|2  = 1 - 


1-1 


1 + 1 


2/3/* 


(8) 


(9) 


(10) 


where  the  value  of  £ is  obtained  from  the  relationships  between  physical  and  auxiliary  coordinates  on  the  interval 
[OB],  see  figure  3: 


and  on  the  interval  f OA\: 


|Z«"  -wX^I) 


d£  , 0 < £ < 1 


(11) 


(12) 


The  value  of  is  the  only  unknown  parameter  left  in  this  problem  to  be  determined  through  the  matching 
procedure. 

It  follows  from  (11)  that  X ~ JV£2/(2V^)  as  X —*  — oo  and  £ — > +oc  and  therefore  velocity  behaves  as  X — > — oo 


dw 

dZ 


V1 


VI 


y/2  0 Xn 


If  one  assumes  that  !+>  = 1,  then  coefficient  K is 

« y/20  Xn  s/2  ft  1 

^ V Vi  7T  v7?  ’ 

see  figure  4 for  JC(ft)  curve,  which  attains  its  maximum  at  3 « 1.734  (or  99.34°). 


(13) 


Figure  4:  Parameter  1C  versus  ft  for  the  simplest  inner  spoiler  problem. 
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2.3.2.  Cavitating  shock  free  leading  edge  Introduce  a local  stretched  coordinate  system  in  the  vicinity  of 
the  wedge-like  leading  edge  Z = X + iY  = z/lrmw.  where  lrmw  -C  a is  the  length  of  the  upper  wetted  portion 
of  the  leading  edge.  As  stagnation  point  coincides  with  the  vertex  of  the  wedge,  value  of  lrmw  is  unknown,  but 
can  be  used  as  a stretching  factor. 

The  flow  pattern  and  corresponding  auxiliary  ((-plane  are  depicted  in  figure  5.  The  flow  region  is  bounded 
by  semi- infinite  line  (OB),  segment  [ OCJ ] and  free  surface  (CD).  Note  that  the  interval  [ OC ] is  of  unit  length 
because  of  the  choice  of  the  stretching  factor.  Like  for  the  inner  spoiler  problem,  the  absolute  value  of  velocity 
at  ‘local’  infinity  and  on  the  free  surface  (CB)  is  an  unknown  parameter  to  be  determined  through  the  matching 
procedure  with  outer  problem. 


Figure  5:  Flow  pattern  and  auxiliary  ((-plane  for  the  inner  leading  edge  problem. 


With  the  correspondence  between  physical  and  auxiliary  planes  (the  first  quadrant  of  the  £ = £ + if/  plane  is 
chosen  as  an  image  of  the  region  occupied  by  the  fluid),  one  can  easily  obtain  the  analytical  inner  nonlinear 
solution  (following  to  Chaplygin  method  of  singular  points)  in  the  form: 


x*(C)  = 


die  - i 

vJJz  = [cTi 

c 


~ih r 


dw 

dC 


N C(C2  + i); 


= >(£if  «■ 


(14) 


The  length  of  the  segment  [OC]  (|0C|  = 1)  is  connected  to  the  unknown  parameter  N through  the  relation 


COl  = ^Jt{1~t2) 

x o 


7/tt 


N 

m = VcQ- 


(15) 


It  is  seen  that  X ~ Nrfl/(4V^0)  as  X — » +oo  and  Y = 0 (£  = 0,  r;  — ► +oo).  Therefore  the  asymptotic  structure 
of  the  inner  solution  at  ‘local’  infinity  as  X — > +oo  and  Y = 0“  is 


die 

d Z 


V 1 


VI 


(16) 


Such  a behaviour  corresponds  to  a well-known  singularity  ‘—1  / 4’  arising  in  the  framework  of  the  linear  cavitation 
theory  at  the  leading  edge  of  the  hydrofoil. 

Note  that  in  the  case  of  7 = 7r  the  inner  problem  reduces  to  a well-known  problem  of  the  flow  in  the  vicinity  of 
the  sharp  leading  edge  of  the  cavitating  hydrofoil  considered  by  Plotkin  [15]  on  the  base  of  a hodograph  method, 
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Figure  6:  Flow  pattern  and  auxiliary  £-plane  for  the  inner  leading  edge  problem  in  the  case  of  7 = 7r. 


see  figure  6.  The  solution  presented  below  is  appeared  to  be  more  efficient  and  convenient: 


, ^=C-i 

A T/2  A *7 


-TT  = N(((2  + 1) 

v^dz  C + 1 dC 


where  L = \OC\  and,  since  X ~ (3/17)  Lrj 4 as  X — > +00  and  77  — + +00, 

(d“'  JP**) 


d Z 


(17) 


(18) 


as  X —*  Tco,  Y = 0 . The  value  of  C and  V//  are  to  be  determined  through  the  matching  procedure. 


2.3.3.  Cavity  closure  region  It  is  obvious  that  outer  linear  solution  loses  its  correctness  in  the  vicinity  of 
the  cavity  trailing  edge  as  z —*  L,  where  L is  the  cavity  extent,  provided  a linearized  closed  cavity  closure 
scheme  (like  Riabouchinsky,  Tulin-Terentev  or  Efros-Kreizel-Gilbarg  models)  is  adopted.  The  outer  expression 
for  the  conjugate  velocity  x°(z)  has  a square-root  singularity  at  2 = L. 

Since  the  latter  linearized  model  is  adopted  in  the  outer  region,  a local  nonlinear  model  is  proposed  of  the 
flow  in  the  close  proximity  of  the  re-entrant  jet,  see  figure  7.  The  upper  half  of  the  flow  pattern  is  shown 
in  stretched  and  rotated  local  coordinates,  the  stretching  factor  being  1 /Se,  where  Se  is  re-entrant  jet  width 
at  infinity  and  rotating  angle  being  je,  where  angle  qe  describes  a direction  of  re-entrant  jet  at  infinity  for 
Efros-Kreisel-Gilbarg  cavity  closure  scheme.  Since  the  flow  pattern  has  axial  symmetry,  the  streamline  (AC) 
is  substituted  by  the  straight  solid  boundary  (ABC).  Then  the  region  occupied  by  the  fluid  is  bounded  by  free 
surface  (AFC)  and  solid  horizontal  line  (ABC),  velocity  vector  at  point  D is  directed  downward.  The  dividing 
streamline  meets  the  solid  wall  at  stagnation  point  B.  This  streamline  subdivides  the  flow  into  two  parts:  the 
first  one  is  the  main  flow  directed  downstream  with  velocity  Yi,  at  local  infinity  as  X — > +00  and  the  second 
one  forms  a re-entrant  jet  of  unit  width  at  infinity  at  point  —00. 

With  the  correspondence  between  the  physical  z = x + iy  plane  and  atixiliary  quadrant  C = C + lrl  shown  in 
figure  7,  the  Chaplygin  method  allows  us  to  write  down  the  following  solution: 


xW 


d w 


C-i 
C + 1’ 


*E  = N£zi. 

dC  c ’ 


(19) 


d Z dw 


dC 


N (C2 

VL  (V  + 2^  + ^-2-5 


where  unknown  parameter  N is  connected  to  the  flow  rate  in  the  jet  at  point  A through  the  relation 


(20) 


N = 


2 VI 


7 r 


(21) 
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Figure  7:  Flow  pattern  and  auxiliary  ("-plane  for  the  inner  cavity  closure  problem. 


It  readily  follows  from  condition  x*(i /)  = i that  / = 1 what  gives  the  distance  L*  between  stagnation  point  B 
and  normal  projection  of  point  F onto  the  horizontal  line  (ABC): 

L*  = jRe  (Z (i))j  = — . (22) 

7T 

Both  Se  and  are  to  be  determined  through  the  matching  procedure  with  outer  linear  solution. 

It  is  seen  that  X ~ (2/7r  as  X — ► +oo  and  ( — > +oo,  that  is  why 


dw 
d Z 


V1 


(23) 


which  corresponds  to  the  square-root  singularity  of  the  conjugate  velocity  for  the  outer  linear  expansion  at  the 
trailing  edge  of  the  cavity  at  z = L. 


2.4.  The  matching  procedure 

It  was  mentioned  above  that  outer  solution  loses  its  correctness  in  the  vicinity  of  the  hydrofoil  and  cavity  edges 
(z  — > 0,  z — » 1 and  z — * L)  where  function  x°(z)  has  singularities  of  1/2  (square  root)  and  1/4  type.  More 
to  the  point,  the  outer  solution  is  not  completed,  for  the  number  of  unknowns  is  greater  then  the  number  of 
conditions  and  coefficient  B in  (1)  or  (2)  is  an  unknown  parameter.  The  main  goal  of  the  matching  procedure 
is  to  close  the  whole  problem  (to  determine  all  the  parameters  in  outer  and  inner  descriptions)  and  to  get  rid 
of  the  singularities  of  the  linear  outer  solution. 

Generally  speaking,  the  matching  procedure  is  carried  out  into  three  steps:  the  first  one  allows  us  to  match 
asymptotic  descriptions  in  the  outer  region  and  in  the  trailing  edge  region  (the  velocity  U\  = in  expression 
(7)  and  the  coefficient  B in  the  equation  (1)  or  (2)  are  determined);  the  second  one  matches  the  outer  description 
with  known  value  of  B and  leading  edge  expansion  to  yield  the  values  of  /w  and  U2  = V//  in  (14)  and  the  third 
one  enables  re-entrant  jet  width  6e  and  velocity  I/3  = V//  in  (19)  to  be  derived.  Notation  Ui,  U2  and  U3  is 
introduced  to  distinguish  velocity  absolute  value  in  the  local  regions. 

First,  it  is  obvious  that 

La3 

( + a 2“  (1  ~x), 

as  x —*  1~,  y = 0“  and  ( — > —a+,  tj  = 0+  (see  figure  2)  and  therefore  the  limiting  form  of  (1)  and  (2)  near  the 
trailing  edge  with  spoiler  is 


Xf(x) 


1 + 


B 

Va(a  + 0 


1 + 


y/2  B 1 
a2^L  'JT^x  ’ 


where  L = cos  2(r/ 2)  > 1 and  a = cot(r/2). 
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Taking  into  account  the  fact  that  X = (x  — l)/e.  where  s = e/l  - is  relative  spoiler  length  this  expression  and 
formula  (13)  are  compared  to  give 

Ui  = Vq  = Vxy/l  + o ; 


„ P T ■ -2  T £ 

B = cos  — sm  — \ — , 

7T  2 2 Mr 


(24) 


where 


de- 


ft is  to  be  underlined  that  all  the  coefficients  A.  B,  C,  T> £k  in  (1)  and  (2)  are  of  the  same  order  of  magnitude. 
That  is  why  s = 0(a2)  as  A = 0(a). 

Second,  it  is  seen  that 

c [* 

^-\L’ 

as  x — i 0+,  y = 0"  and  £ — 1 ► 0~.  therefore  the  limiting  form  of  the  classical  solution  (1)  near  the  leading  edge  is 

xf(x)  ~ 1 + h/f  (-  / °(f)  At  + - +c)  ~ 1 + i/4  (~  f m dt  + B + oc)  4r  > 

V C U J Jt{a  - 1)  a ) V a2  1 7T  ,/  Jt(a  - 1)  ) 


and  for  nonquadrature  solution  (2)  is 


a (B 


X2l(x)  ~l+iJ-^-+C-2(S1+T>2-S3-Vi+£6+Ve-.. 

~ 1 + ^ jj.  {B  + aC  - 2 a {£,  + V2  - £,  - VA  + £s  + V6  -...)}  , 
the  value  of  the  coefficient  B being  given  by  (24). 

The  matching  procedure  for  these  expansions  and  expression  (16)  considering  the  fact  that  Z = z/(lw/l)  results 
in 


L ir4Q  L 

J 4^4  q2 


where 


U2  = Vo  = VxVT+a  ; 

+ <xC  — 2 a {£\  + — £3  — TT  + £5  + 2?6  — • • •)  J 

y/ir 


(25) 


Note  that  T>k  = £k  = 0,  k = 1, 2, . . . in  the  case  of  the  supercavitating  shock  free  flat  plate.  Expressions  (25) 
show  the  influence  of  flow  parameters  including  spoiler  geometry  onto  the  length  of  the  upper  wetted  portion 
of  the  leading  edge  lw  for  a given  value  of  angle  7.  It  follows  from  the  latter  equation  (25)  that  lw/l  = O (a4) 
as  C.  T>k,  £k  are  of  the  order  of  0(a)  and  B = Q(\/g)  = O(a). 

Third,  in  the  vicinity  of  the  cavity  trailing  edge  x — * L~ , y = 0+  and  £ — > +00,  r)  = 0+  (see  figure  2) 


L — x 

wherefrom  the  limiting  form  of  outer  solution  (both  classical  and  nonquadrature)  is 


X31  (x)  ~ 1 + i.A|  ~ 1 + L4 


L — x 


Taking  into  account  the  fact  that  X = (x  — L)/(Se/l),  where  8e/l  is  a half  of  the  relative  width  of  the  re-entrant 
jet,  this  expression  along  with  (23)  are  compared  to  yield 


c 

U3  = v0  = Too/IT^  ; (f  = ^LA2. 


(26) 
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The  expressions  show  relation  between  the  width  of  the  re-entrant  jet  and  flow  parameters  (including  spoiler 
length  e and  inclination  angle  /?).  Note  that  Se/l  = 0(a2)  because  of  A = 0(a). 

The  ‘information  stream’  in  the  matching  procedure  is  directed  from  the  trailing  edge  of  the  hydrofoil  to  the 
leading  one  and  to  the  trailing  edge  of  the  cavity  as  well.  In  fact,  the  spoiler  geometry  parameters  0 and  e 
dictate  the  value  of  the  coefficient  B which,  in  turn,  defines  the  leading  edge  and  re-entrant  jet  characteristics. 
The  additive  composite  solution  for  the  conjugate  velocity  can  be  constructed  in  a following  manner: 

Xc(z)  = X°(z)  + x! (z)  + xtiz)  + Xjs(z)  - X? (z)  ~ xf(z)  - xf(z)  , (27) 

where  subscripts  1,  2,  3 correspond  to  local  flow  problems  in  the  vicinity  of  the  spoiler,  leading  edge  and  cavity 
trailing  edge.  Substituting  composite  solution  xc  instead  of  outer  x°  into  formulae  (3)-(6)  enables  one  to  derive 
hydrodynamic  coefficients  which  are  everywhere  valid  and  to  calculate  correct  flow  pattern. 


2.5.  Exact  solution  to  the  nonlinear  flat  plate  problem  and  asymptotic  analysis 


Consider  a nonlinear  problem  of  the  flow  around  a supercavitating  shock  free  flat  plate.  It  is  useful  as  a 
verification  of  the  asymptotic  results  obtained  in  the  previous  sections.  The  problem  under  consideration  is 
that  of  the  theory  of  jets  in  an  ideal  fluid  and  has  to  be  treated  by  corresponding  methods.  The  physical  2 = x+vy 
and  auxiliary  C planes  are  shown  in  figure  8.  The  same  notation  is  used  as  in  section  2.1  where  the  general 
shock  free  cavitating  problem  is  formulated.  The  velocity  absolute  value  on  the  free  surfaces  is  vq  = Voo  \/l  + a . 
where  a > 0 is  cavitation  number.  The  region  occupied  by  the  fluid  Ls  bounded  by  solid  segments  [AC],  CO] 
and  [OB]  and  cavity  surfaces  (BE)  and  ( AE ),  point  E is  that  at  infinity  (second  Riemann  sheet).  Point  at 
infinity  F (first  Riemann  sheet)  has  the  image  (jo  = b + ic  and  stagnation  point  F - image  Coo  = d + if  on  the 
first  quadrant  of  the  auxiliary  C-plane. 


Figure  8:  Flow  pattern  and  auxiliary  C-plane  for  the  nonlinear  shock  free  flat  plate  problem. 

With  the  correspondence  between  the  planes,  the  Chaplygin  method  allows  us  to  write  down  the  following  exact 
solution: 


xn(C)  = 


d w 
Vodz 


C ~a 

C T u 


~th 


C-A^C-CoHC-Co) 


C + tj  (C  + Co)  (C  + Co)  ’ 


dw  ,TC(C2-«2)(C2-C02)(C2-C(!) 
dc  (C2  + i)(C2  - C2o)2CC2  - CSo)2  ’ 


wherefrom 


(28) 

(29) 


< 


z(f)  = ^ei(7-)  [ £(C  + q)  /7r(C~«) 

Vo 


1+7/7r  (f  U + t 

(C2  + 1)  VC^i 


(C  + Co)2(C  + Co)2  u 

(c2-ej2(c2-cio)2 


(30) 
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Five  conditions 

zc  = z(a)  = l el(?r~“) ; zb  = z(0)  = ; za  — zc  = z( oo)  — z(a ) = lw  el(7~“^ ; 

jf£d{  = o 

Wo  / dC 

<oo 

generate  a system  of  seven  nonlinear  equations  in  eight  real  unknowns  a,  t,  (o  = b + ic.  Coo  = d + if,  N and  lw. 
An  additional  eighth  condition  is  to  be  used  to  close  the  problem,  namely  connected  with  direction  of  the 
re-entrant  jet  at  point  E.  Assume  that  angle  /i  is  given,  then  the  condition  is 

Xn(i)  = e~iM, 


0(1)  = IE 


(31) 


where  9{rj)  denotes  direction  of  velocity  vector 

2 7 a b b 2/3  t 

Oyn)  = — argx(w?)  = — a + 7 arctan 2 arctan 2 arctan arctan  — 

7T  77  r)  — c T]  + C 7T  T) 

The  pressure  distribution  coefficient  is 

Cp(£)  = 1 - ixn(C)l2  = 1 


c - a 

27/yC2-26£  + 62+c2\2 

z-t 

20/n 

£ + a 

U2  + 25£  + b2  + c2  ) 

z+t 

(32) 


(33) 


where  £ = Re  (C)  is  connected  to  z through  relation  (30).  Lift  and  drag  coefficients  are 


CF  = CD+  i CL 


-j(l+*) 


/ 


d w d w , , 

dcd2‘i{ 


(34) 


On  the  other  hand,  using  residue  theory,  one  arrives  at  the  following  relationships  (which  are  correct  for  an 
arbitrary  cavitating  hydrofoil): 


n 2 q (,  v 0 
C d = : — j I 1 - - — cos  fi 

^OO  ^ \ ^OO 


2 Q 


Vo 


CL  = — 7 — sin  /j, 

Woo  I \q  Woo 


(35) 


where  q = 2 ScVq  denotes  the  flow  rate  in  the  re-entrant  jet  E (Se  is  a half  of  the  width  of  the  jet)  and  F is 
circulation  along  a large  contour  completely  surrounding  the  cavitating  foil  and  the  cavity  and  enclosing  most 
of  the  flow.  Note  that 

/ (IF  ^ ^ 

res  = J)  — dC  = i + iq 


and.  moreover, 


q = 7 tN 


(1  + a2)  ((1  + b2  - c2)2  + 462c2) 


2 ((1  + d2  - f 2)2  + 4 d2PY 
Another  form  of  the  force  coefficient  CF  = Co  + iCi,  is 


Cf  = ——  {—  v^l  + (T  e‘MRe  (res)  + res)  . 

Woo  \ / 


(36) 


Let  us  analyse  the  solution  to  the  problem  under  consideration  in  the  case  of  a — ► 0,  a — *■  0,  e/l  — > 0 
and  lw/l  —>  0,  that  is  under  assumption  that  hydrofoil  brings  small  perturbations  into  the  inflow.  Thorough 
asymptotic  analysis  of  the  problem  is  given  in  [1]  and  below  just  final  results  are  shown. 

It  is  clear  from  what  was  done  above  that  a — * +00,  Co  “ 1 * i (b  —> • 0 and  c — » 1)  and  f — > 0 if  /?  = 0(1)  and 

7 = 0(1). 
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Taking  account  of  only  linear  terms  in  equations  (28)-(32)  under  those  circumstances  results  in  (c  — 1)  <C  b (if 
/.*  = 7 r + 0(a))  and 

. . 8/Mo  1 fe  L-n 

a ~ 4jdoa  + - 


Q 


l 4 74  (/2  +d2)2 


7T  1 


16  /2d2  \/  2 doa  - fd°  ^ 


dw 

v~>dz 


l-^-46 


C W t 

C2  + 1 7T  ( ’ 


do  — + /«  + 


W 

7 r 

dw 


P f /£ 

7T  -/R  V l 

1 


(37) 


v/R  V / 

—Na2C(£2  + 1) 


dC  (C4  + 2C  2 + (/2+dg)2)2' 


(38) 


where  Ln  denotes  ‘nonlinear’  cavity  extent,  z = exaz  and 

l 

d0=  ' 


7 /it 


d£. 


o o 

The  conformal  mapping  of  the  first  quadrant  of  £-plane  onto  the  upper  semi-plane  w has  the  following  form 

where  uQ  = 2dof  ~ Ln/l . Substituting  this  expression  into  (38)  gives 
dw  „ . 2/3  too 


(39) 


1 - i 


7T  y'u  (u.  + -U0) 


+ i 


2u 


4 b 27  1 2(5  \ 

— u 1 1 

(40) 

uq  tt  a 7 r y 

(41) 

<lz  ,r  „ 

Bearing  in  mind  that  ‘linear’  L and  ‘nonlinear’  L„  cavity  extent  are  connected  to  each  other  through  the  relation 
(due  to  a shift  of  coordinate  systems) 

l = 1 + t 

and  L = cos~2(t/2),  we  find  that 


sin  - 


7 r — r 


cos  - 


7 r — r 


do 


T 

Sill  — 
2 


■;  /■ 


Substituting  these  limits  into  asymptotic  expansions  obtained  above  gives  expressions  coinciding  with  those 
reduced  from  outer  (2)  and  composite  (27)  solutions  in  the  case  of  a supercavitating  shock  free  flat  plate  when 
f(x)  = — ax  and  T>k  = £k  = 0.  Moreover,  the  asymptotic  and  nonlinear  approaches  give  the  similar  formulae 
for  the  width  of  re-entrant  jet  28e  and  upper  wetted  portion  of  the  leading  edge  lw/l  (compare  (26)  and  (25) 
with  (37)). 


2.6.  Numerical  results 

Mathematics  for  Windows  computer  mathematical  environment  was  used  to  obtain  all  the  numerical  results 
shown  in  this  section. 

Flow  pattern  and  hydrodynamic  coefficients  are  given  in  figure  9 for  nonlinear  theory  (NLT)  and  matched 
asymptotics  (MAE)  for  a shock  free  flat  plate.  The  flow  parameters  are  chosen  to  be  a = 10°,  a = 0.5,  e = 0.05, 
P = 90°,  7 = 60°  and  /i  = 200°.  The  position  of  stagnation  point  F is  also  shown. 

Pressure  distribution  coefficient  Cv  is  depicted  for  such  a flat  plate  in  figure  10  (for  the  same  set  of  the  flow 
parameters).  Numerical  results  for  the  linear  theory  (outer  asymptotic  expansion  x°(z))  is  shown  as  well. 


